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Ocean rogue waves are large and suddenly appearing surface gravity waves [1], which may cause 
severe damage to ships and other maritime structures [HIS]. Despite years of research, the exact 
origin of rogue waves is still disputed [4]. Linear interference of waves with random phase has often 
been cited as one possible explanation mi. but apparently does not satisfactorily explain the 
probability of extreme events in the ocean [Till]. Other explanations therefore suggested a decisive 
role of a nonlinearity in the system [2113 EHII]- Here we show that linear interference |12[ 113) of a 
finite and variable number of waves [nnn] may very well explain the heavy tail in the wave height 
distribution. Our model can explain all prototypical ocean rogue waves reported so far, including 
the “three sisters” m as well as rogue holes We further suggest nonlinear time series analysis 
|18II19| for estimation of the characteristic number of interfering waves for a given sea state. If ocean 
dynamics is ruled by interference of less than ten waves, rogue waves cannot appear as a matter 
of principle. In contrast, for larger numbers, their appearance is much more likely than predicted 
by parameter less models [3 or longterm observation [ 8 ]. The pronounced threshold behavior of our 
model enables effective forecasting of extreme ocean waves. 


Surface gravity waves are a result of wind blowing over 
the ocean surface [T] . As wind directions differ across the 
ocean and in time, waves with different origin may start 
to interfere with each other at a given fixed location. In 
narrow confined waters like small ponds, one can some¬ 
times observe the build-up of perfectly plane wavefronts 
across the entire surface width. In such a case, virtu¬ 
ally no interference appears. In the much larger ocean 
system, in contrast, one anticipates an inextricable inter¬ 
ference of many waves, which is typically treated in the 
limiting case of an infinite number of interfering waves 
011]. Mathematically modeling this interference process 
indicates a heavytail distribution of wave heights. Con¬ 
sidering only linear interference, a Rayleigh distribution 
results |1|. Accounting for nonlinearities of the ocean 
system, the Tayfun distribution with an even more pro¬ 
nounced heavy tail emerges m- Unfortunately, both 
distributions tend to significantly overestimate the ap¬ 
pearance frequency of rogue events in the ocean IZIIHI- In 
the following we will show that the heavy tail of the distri¬ 
bution can be adjusted by assumption of a finite number 
N of interfering waves [III [T5|, as illustrated in Fig. 

In particular, if N drops below a certain threshold value, 
rogue waves become completely impossible whereas their 
likelihood rapidly increases with growing N above the 


threshold. Moreover, the number N is shown to be vari¬ 
able, depending strongly on location and weather condi¬ 
tions. Consequently, estimation of N enables early warn¬ 
ing of situations that are prone to rogue wave formation. 

Nonlinear time series analysis offers a method that al¬ 
lows effective estimation of the number of interfering el¬ 
ementary waves from observations at a single location 
in the ocean. One suitable method is the Grassberger- 
Procaccia analysis (GPA) [THl US] , and it delivers an es¬ 
timate for the dimension D of the parameter space of 
the system. While the dynamics of simple mechanical 
oscillators is fully characterized by knowledge of a single 
phase, the interference of many oscillators at one point 
on the surface of the ocean, in contrast, is ruled a large 
number D = N oi individual phases, see Figs. g,h,c). 
Previous studies of ocean wave dynamics indicated val¬ 
ues of D reaching from 7 to 10.5, with a certain tendency 
of lower values appearing in calmer waters HmH]. Ap¬ 
plying the GPA to wave height measurements recorded 
on January 1, 1995 on the Draupner platform [23] . we 
determine a dimension D in the range of 12 to 13, see 
Fig. [^a). The underlying measurement [Fig. |^b)] is the 
first record of a rogue wave, which ultimately confirmed 
their existence. The significant wave height Hs in this 
event was 12 m, and the rogue wave exhibited a wave 
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FIG. 1: Simulation of the interference on a surface with a normally distributed angular spread of N waves amounting to 
SOmrad. In this simulations, a nonlinear correction was included to reduce the depth of the troughs and steepen the wave 
crests, (a) Assuming random phases between the waves. A highly random wave pattern emerges with no wave significantly 
exceeding the significant wave height. N = 20. (b) Simulation of a rogue hole, with synchronized phases adding up to negative 
interference at the central point of the propagation axis. N = 40. (c) Same with phase flipped by tt and N = 20, resulting in 
a single wave with extreme height > 2.5iFs Depending on the diffusion coefficient, leading and trailing troughs appear, similar 
to the observation of the Draupner event [23] . (d) Same with weaker phase diffusion. A rogue wave with structure reminiscent 
of the “three sisters” m emerges. 


height of 25.6 m trough to crest. For comparison, we also 
ran the GPA on other data sets [Mj recorded during the 
Draupner storm, which provided a nearly identical result 
for D. It should be noted that significant wave heights 
were considerably smaller (Hg < 4 m) in previous stud¬ 
ies mnn. The clearly deterministic signature of ocean 
wave records [21] suggests that the dynamics of the wa¬ 
ter surface arises due to interference of a relatively small 
number of elementary waves. Moreover, this number is 
certainly not a constant and may also depend on factors 
other than Hg- 

Linear interference of waves with random phase is a 
well-known assumption for modeling ocean waves [6]. 
Similar statistics are relevant for the amplitude distri¬ 
bution of multiply scattered radar signals [H] . Recently, 
the emergence of rogue wave statistics due to purely lin¬ 
ear interference of optical waves was observed in two dif¬ 
ferent experimental situations [Hdg, and similar linear 
rogue waves were seen in microwave experiments m- 
The experiments by Liu et al. beautifully confirmed the 
Rayleigh distribution suggested for the limit of an infinite 
number of interfering waves [25] . In all these examples, 
the fundamental physics of the random interference phe¬ 
nomenon appears nearly identical, yet with a differing 
number N of interfering waves. Figurej^a) shows simula¬ 
tions of randomly interfering waves on a two-dimensional 
surface in comparison to the situation of local perfect 
constructive interference of the waves, see Figs, l(b-d) 
and movie sequence in the SI. Depending on the rate of 
the modeled phase diffusion process [25] , either the char¬ 
acteristic preceding and trailing troughs of the Draupner 


event emerge [Fig. 1(c)], or a pattern appears that is 
reminiscent of the famous “three sisters” rogue wave [T5] 
[Fig. 1(d)]. Moreover, model parameters can also be ad¬ 
justed to produce an isolated rogue hole [T7j [Fig. 1(b)], 
which appears to be the most mysterious rogue wave vari¬ 
ant reported to date. 

Figure]^ shows results of more detailed numerical sim¬ 
ulations of the random interference of a finite number N 
of waves, see Methods sectio n for details. The result¬ 
ing mean height amounts to '/ttN /2, see Fig. [^a). The 
resulting probability densities pN{h) are skewed, and a 
heavy tail starts to form with increasing N. We adopt 
the fairly strict definition that requires rogue waves to 
exceed the significant wave height Hg by a factor 2.2 and 
find that this condition can only be met for N > 10, see 
Fig. [^a). Using TV = D = 12 as suggested by the GPA, 
the probability for rogue wave appearance amounts to 
P(/i > 2.2Hs) « 2 X 10“®, which agrees favorably with 
an exhaustive analysis of 10® days of wave height obser¬ 
vations [H], see Fig. [^b). For example, assuming a wave 
period of 15 s, these estimates translate into the appear¬ 
ance of a single rogue event within hundreds of days of 
permanent Draupner storm conditions. Going to signifi¬ 
cantly higher values oi N ^ 20, the probability for rogue 
wave dramatically increases, with resulting probability 
density functions slowly converging towards a Rayleigh 
function. Under these severe storm conditions, probabil¬ 
ities for the emergence of rogue waves may readily be ten 
times higher than during the Draupner event (iV = 12). 
Comparing the simulated probabilities with the analy¬ 
sis of Christou and Ewans [S], nevertheless, it appears 
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FIG. 2: (a) Dimensional analysis of the Draupner wave record 
with the Grassberger-Procaccia method (see Methods sec¬ 
tion). Shown is the correlation sum slope as a function of 
embedding dimension for the original Draupner record (green) 
and a second measurement from the same campaign (blue). 
Both records contain a rogue event and exhibit a plateau of 
the slope for embedding dimensions ranging from 70 to 130. 
This results in an estimate for the information dimension D of 
the record of 12 to 13, see pink shaded area, (b) Reconstruc¬ 
tion of the Draupner event using the model assumption of 
an interference of 12 waves of equal and constant amplitudes 
but temporally varying phases <Pj(t), {j = 1, ■■■, 12), (original 
data: black curve, symbols: reconstruction), (c) Phase func¬ 
tions ifijit) determined by simulated annealing. Successful 
reconstruction indicates that a completely linear explanation 
of rogue wave formation via phase diffusion of a small number 
of elementary waves is possible. 

striking that the observed trend of wave heights does 
not follow the simulated probability for any fixed number 
N. This finding reflects the fact that this extensive wave 
record was not recorded under constant weather condi¬ 
tions and therefore corresponds to a mixture of probabil¬ 
ity distributions with various values of N. In fact, the 
comparison suggests that the ocean dynamics is proba¬ 
bly well characterized by TV = 6 throughout most of the 
time. However, occasional storms seem to even exceed 


N = 25. A second explanation for the largest observed 
rogue events h > 2.5Hs in the Christou record may be 
nonlinear steepening or compression effects, which fur¬ 
ther amplify the heavy tail, similar as in the Tayfun cor¬ 
rection [50] to linear statistics. Such mechanisms have 
been previously discussed as super rogue waves m- 

The variable character of N appears key to under¬ 
standing the rogue wave phenomenon. As the waves are 
generated by winds across the ocean surface, one expects 
to see an increase of the average number of interfering 
waves, e.g., when wind directions change over a wide 
range or when turbulence sets in. A strong variability 
of wind directions is expected to generate waves within a 
wider angular range and may therefore give rise to larger 
values of N. Moreover, an increasing build-up of large 
waves on the ocean surface will also serve to increase the 
turbulent nature of the wind blowing across this undu¬ 
lated surface. The hypothesis of a connection between 
strong winds and frequent rogue waves is also well con¬ 
firmed by the observations of Christou and Ewans [S] 
who found hundreds of rogue waves at strong winds, but 
could not find a single one at below-average wind con¬ 
ditions. Additionally, they also observed an increase of 
rogue events for the case of a wide angular spread of the 
swell. Let us further remark that the number of inter¬ 
fering waves may also be influenced by particular coastal 
shapes, subsea topography, or ocean currents, which may 
have a focusing or cumulating effect on otherwise spa¬ 
tially separated waves [2S1HH] ■ In order to identify such 
possibly dangerous regions, we suggest using the correla¬ 
tion dimension D for characterization of the sea state. As 
indicated by the example analysis of the Draupner data, 
relatively short time series of 20 minutes length already 
suffice for an estimate of D, and required computational 
times are on the order of a few minutes. Values above 
the threshold of I? « 10 indicate the possibility of rogue 
waves, with a rapidly increasing probability for higher D. 
At open waters far away from any coastline, it may addi¬ 
tionally prove useful to carefully monitor the wind situ¬ 
ation, in particular the variability of the wind directions 
or the angular spread of the swell. A further possibil¬ 
ity for identifying rogue wave situations may be given by 
measuring the transverse extent of the wave crests, which 
is expected to shrink with increasing angular spread of 
interfering waves. 

In conclusion, our analysis indicates that rogue waves 
in weakly nonlinear systems may simply appear due to 
linear interference of a finite number of plane waves larger 
than ten. No nonlinearity in the oceanic system is there¬ 
fore required to explain the formation of a heavy tail, yet 
nonlinearities will certainly give rise to modification of 
the probability density functions, and they may also ex¬ 
plain the extremely steep waveforms reported for oceanic 
rogue waves. Our findings indicate particular wind con¬ 
ditions as the real origin of rogue waves in the sea, and 
nonlinearities in the atmospheric system certainly play a 
decisive role in generating strongly varying wind direc¬ 
tions. Quite surprisingly, moderate wind conditions may 
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FIG. 3: Simulation of statistical distributions resulting from 
the interference of N waves with unity amplitude and random 
phases. For each value of N, 10® resulting wave heights have 
been computed, i.e., about 4 times more than available in 
Ref. [8]. (a) Resulting average wave height, significant wave 
height Hs as well as maximum observed wave height as a 
function of N. For identical phases, a maximum possible 
wave height results at N times the height of the individual 
input waves. Using the threshold condition of 2.0 or 2.2Hs, 
no rogue waves are possible below N = 8 or 10, respectively, 
(b) Computed probabilities of exceeding a given wave height 
in units of Hs- Cases from = 6 to 100 are shown as lines. 
Observations from Christou and Ewans are shown as dots 
[^. The Foristall distribution is shown as a dashed line [3. 
The fact that no single simulation fits to the observations is 
indicative of a strong variation of N due to changing weather 
conditions during the observations. 


completely thwart the formation of rogue waves whereas 
the fully developed turbulence in gale-force storms may 
invoke these “monsters of the deep” at much higher prob¬ 
abilities than previously thought possible. The system 
size estimate D from the GPA analysis may play a piv¬ 
otal role in forecasting rogue wave situations, and it may 
further help to identify particularly dangerous areas in 
the ocean. Finally, given the critical dependence on wind 
conditions, even very minor climate changes |30) strongly 
affect the probability of rogue waves. A 10% increase of 
wind speeds in an average storm may readily double the 
appearance of rogue waves. As the bottom line, the an¬ 
swer to the ocean rogue wave mystery does not seem to 
lie in the depth of the ocean as is commonly believed; in¬ 


stead, this answer appears to be literally blowing in the 
wind. 


Methods 


Probability Density Functions 


Various probability distributions been discussed for 
ocean waves. In the case of an infinite number of lin¬ 
early interfering waves, a Rayleigh density function 

4'Rayleigh(^) = ^ 

of wave heights h results [6] . cr is a scale parameter, which 
parameterizes both the average height and the standard 
deviation of the waves. Taking nonlinearities of the ocean 
wave system into account, M. A. Tayfun derived the fol¬ 
lowing correction |20] to the Rayleigh distribution 


TTayfun(^) — 2/3 I 1 


a/1 + ‘2h/B 


X exp 


-r Vl + 



( 2 ) 


where 


/? 


1 


(3) 


and 7 is a nonlinearity parameter, which is typically cho¬ 
sen such that the product ycr is in the range from 0.1 
to 0.3, further increasing the heavy tail in the distribu¬ 
tion. While both of these distributions have been re¬ 
ported to overestimate the heavy tail, best agreement 
with observed data was reported for 


PForristall(/l) = ^ exp 


( /i2.126 X 

2.105 0-2 j’ 


(4) 


as was empirically found by G. Z. Forristall [S]. The ex¬ 
cellent agreement with the total of 30 years of wave data 
was independently confirmed by Christou and Ewans 
[5]. Using any of the above definitions, the probability 
V{H > h) for the wave height H to exceed a given value 
h is computed as 


P(iJ > h) = J p{s)ds j J p{s)ds. 


(5) 


Grassberger-Procaccia Dimensional Analysis 

Neglecting dynamical changes in the amplitudes, the 
dynamics of the ocean system at one fixed location is 
ruled by a number of N individual phases, which is 
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identical to the dimension D of the parameter space, 
i.e., D = N. The number D can be estimated by the 
Grassberger-Procaccia algorithm 1 18) (GPA), which is a 
standard method of nonlinear time series analysis. The 
GPA analyzes a time series of a physical quantity x. This 
time series has a length N and is discretely sampled with 
a constant sampling rate. We denote the time series as 
X = {xi,X 2 ,x^,... ,xn}- From this time series x, sub¬ 
series iji = {xi, Xi+i, Xi+ 2 , ■ • ■, Xi+rn\ of length m are se¬ 
lected. The length m is called the embedding dimension. 
Each sub-series yi is compared to all sub-series yj with 
j > i hy calculating Euclidian distances 


njm =11 Vz - Vj 11 = 


2+771 


.^\xk-Xk+j-^\^■ ( 6 ) 

\ k=i 


For m ^ N the Euclidian distances Vijm are accumulated 
in the correlation sum 


Cm{r) = - njm), 




(7) 


j>i 


where 9{r) is the Heaviside step function. In order to 
avoid false interpretation due to detection noise and lim¬ 
ited sample length, Cm {x) is only analyzed at the interval 
[0.01 T’maxjO.l Tniax]) where Tmax Is the largest Euclidian 
distance found in the data series. The correlation sum 
Cm{x) increases monotonically according to 

Cm{r) oc (8) 


with the exponent v. For small values of m, v increases 
with embedding dimension m until it reaches a satura¬ 
tion value dsaf For a range of values of m, a plateau 
with nearly constant dm = dsat is typically observed 
(cf. Fig. [^, provided a sufficient number of data points 
in the analysis. Strictly speaking, the correlation dimen¬ 
sion dsat delivers a lower bound for the (Hausdorff) di¬ 
mension D of the parameter space of the svstem[TB]. In 
many cases, dgat directly delivers a good estimate for the 
dimension D. In case of a number of N interfering me¬ 
chanical oscillators of constant amplitude, the parameter 
space of the system is composed of N independent phases 
LpN- In this situation, we can set D — N > dsat- 


Superposition of Waves as a Two-Dimensional 
Random Walk Problem 


The problem of the resulting amplitude of N interfer¬ 
ing waves is mathematically equivalent to computing the 
probability of the length h of two-dimensional random 
walks 0 HU [TS] involving N steps 


h = 


N 

^/ife exp(j(/jfc) , 
k=l 


(9) 


where tpk € [0, 27r[ are univariate uncorrelated random 
phases and hk are the heights of the individual waves. As¬ 
suming hk = hQ, the absolute maximum possible result¬ 
ing wave height is iV/ig, and for sufficiently large N the 
mean height is « y/'xNho/2, cf. Fig.^a). The threshold 
for rogue waves scales as « 3.1/io-\/fV7lf all ipk are identi¬ 
cal, the maximum possible amplitude of hma.x = Nho oc¬ 
curs. Gomparing the rogue threshold with the maximum 
possible wave height, it becomes immediately clear that 
a minimum of 10 linearly interfering waves is required 
to exceed the rogue threshold. The resulting probability 
densities pnih) are skewed, and a heavy tail starts to 
form with increasing TV. Computing h for a large num¬ 
ber of randomly chosen pk enables the computation of 
probability density functions for a given TV. For the ex¬ 
amples shown in Fig. a total of 10® random sets has 
been employed for any of the shown values of TV. 


Stochastic Model of Rogue Wave Emergence 

In numerical simulations, we assumed linear interfer¬ 
ence of TV elementary waves [6] of unity amplitude accord¬ 
ing to 


N 

h{t, r) = '^2 cos (^kjr + cot + (pj(x,t)"j . (10) 

T=i 

Here a constant angular frequency of the waves oj has 
been assumed, t is the time, and r = (x, y) is the 
spatial coordinate. The wave vectors are defined as 
kj = {cos{0j),sin{9j)), with a Gaussian distribution of 
random angles 9j of vanishing mean. The standard de¬ 
viation A9 of the angular distribution of waves deter¬ 
mines the lateral extent of a rogue wave along the y di¬ 
rection. Using a 27r univariate angular distribution in¬ 
stead, one can simulate the emergence of circular sym¬ 
metric rogue waves as observed by Arecchi et aim- A 
one-dimensional diffusion process of the phase functions 
ipj{x,t) along the propagation direction[2ni is computed 
at equidistant positions Xi = iA^ according to the recur¬ 
sive law 


iPj{xi±i,t) = (pj{xi,t) -y ( 11 ) 

with a diffusion coefficient and a set of independent 
Gaussian random phases tpij with mean zero. This model 
assumption induces a random phase walk, which never¬ 
theless leaves the average frequency uj unaffected. The 
characteristic length or time over which dephasing ap¬ 
pears is related to the prediction length or linear correla¬ 
tion time [24]. Rogue waves appear by setting ipj = 0 at 
t = 0 and X = 0; rogue holes can be generated by forcing 
all ipj = TT. 

The above model can also be inverted to find a set of 
phase functions Pj{t) that give rise to a measured wave 
record h{t), cf. Figs. ib, c). Starting at the time t = 0 
of the rogue event, we used a simple simulated annealing 
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variant to adjust the individual phases ^j(t) in a round- 
robin fashion to the end of obtaining best agreement with 
the measured record of the Draupner event. This proce¬ 
dure is repeated for each time step in the measured time 
series. As this phase retrieval is an ill-posed problem, 
rapid temporal oscillations of the phase functions may lo¬ 
cally appear. These oscillations can be mostly suppressed 
by a suitable penalty term, with the noted exception of 
the deep trough areas of the Draupner wave. This may 
be an indication for unaccounted nonlinear shaping in the 
immediate vicinity of the rogue wave. It should also be 
noted that phase retrieval is certainly highly ambiguous 
in this situation, but nevertheless indicates that rogue 
wave formation can be explained due to phase diffusion 
in the interference of a hnite number of elementary wave 
with constant and equal amplitude. The phase dynamics 
near t = 0 shows essentially a linear dephasing behavior 
with time, mostly being caused by differing frequencies of 
the interfering waves. For the Draupner event, the phase 
diffusion rate between the individual waves amounts to 
47r per minute. 

Temporal phase diffusion can be simulated in a com¬ 
pletely analogous way {U = iAt) by 

(f^) 

with Gaussian random phases (pij and an independent 
diffusion constant Dt- In Fig. N has been set to 20, 
which exceeds the estimate from the dimensional analysis 
{D = 12 to 13). For simulation of a Draupner-like event, 
the diffusion constant was chosen such that dephasing 
appears on the scale of one oscillation period according to 
the above consideration. For the “three sisters,” smaller 
values of need to be used, giving rise to dephasing af¬ 
ter more than one oscillation cycle. Provided a sufficient 
number of rogue wave data sets, the exact value of the 
diffusion constant can be estimated from a linear corre¬ 
lation analysis of the wave record. The angular spread 
of the 9j has been adjusted to 50 mrad to allow a certain 


vertical extension of the rogue waves. Narrowing down 
the spread, wider and wider walls of water appear. Pa¬ 
rameter Dt, finally, determines how long a rogue wave 
will survive while propagating on the ocean surface. At 
the current stage, we are not aware of any measurements 
that would allow estimating this parameter. 

Accounting for the nonlinearity of the ocean system, 
one can include a nonlinear correction in the simulation, 
which effectively reduces the depth of the troughs and in¬ 
creases the crests. When doing so, a substantially larger 
number N of interfering waves is required to form a rogue 
hole whereas the likelihood of rogue holes and waves is 
equal in a completely linear model without correction. 

Our model relates the emergence of ocean rogue waves 
to only three parameters, namely N, D^p, and Dt- The 
effective number N of interfering waves can be estimated 
via the GPA analysis, the phase diffusion coefficient 
can be related to the characteristic number of crests and 
troughs in a rogue event or to a linear correlation analysis 
of the time series. The phase diffusion rate in time is cur¬ 
rently unknown, but this process is expected to be slower 
than the primary phase diffusion, as rogue waves would 
otherwise only remain stable for a few seconds. Among 
the three mentioned parameters, N is the most important 
one, since its knowledge allows predicting probabilities of 
wave heights by isolating the dynamics at one fixed point 
on the ocean surface. 
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